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We analyze an ideal gas like models of a trading market. We propose a new 
fit for the money distribution in the fixed or uniform saving market. For the 
market with quenched random saving factors for its agents we show that the 
steady state income (m) distribution P(m) in the model has a power law tail 
with Pareto index v exactly equal to unity, confirming the earlier numerical 
studies on this model. We analyze the distribution of mutual money difference 
and also develop a master equation for the time development of P(m). Precise 
solutions are then obtained in some special cases. 

1 Introduction 

The distribution of wealth among individuals in an economy has been an 
important area of research in economics, for more than a hundred years. Pareto 
[1] first quantified the high-end of the income distribution in a society and 
found it to follow a power-law P{m) ~ m _ ( 1+I/ ' , where P gives the normalized 
number of people with income m, and the exponent v, called the Pareto index, 
was found to have a value between 1 and 3. 

Considerable investigations with real data during the last ten years re- 
vealed that the tail of the income distribution indeed follows the above men- 
tioned behavior and the value of the Pareto index v is generally seen to vary 
between 1 and 2.5 [2, 3, 4, 5]. It is also known that typically less than 10% of 
the population in any country possesses about 40% of the total wealth of that 
country and they follow the above law. The rest of the low income popula- 
tion, in fact the majority (90% or more), follow a different distribution which 
is debated to be either Gibbs [3, 6] or log- normal [4]. 

Much work has been done recently on models of markets, where economic 
(trading) activity is analogous to some scattering process [6, 7, 8, 9, 10, 11, 
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12]. Wc put our attention to models where introducing a saving factor for 
the agents, a wealth distribution similar to that in the real economy can be 
obtained [7, 8] . Savings do play an important role in determining the nature of 
the wealth distribution in an economy and this has already been observed in 
some recent investigations [13]. Two variants of the model have been of recent 
interest; namely, where the agents have the same fixed saving factor [7], and 
where the agents have a quenched random distribution of saving factors [8] . 
While the former has been understood to a certain extent (see e.g, [14, 15]), 
and argued to resemble a gamma distribution [15], attempts to analyze the 
latter model are still incomplete (see however, [16]). Further numerical studies 
[17] of time correlations in the model seem to indicate even more intriguing 
features of the model. In this article, we intend to study both the market 
models with savings, analyzing the money difference in the models. 

2 The model 

The market consists of N (fixed) agents, each having money m,i(t) at time 
t (i = 1,2,..., N). The total money M (— J^. m,i(t)) in the market is also 
fixed. Each agent i has a saving factor A; (0 < A, < 1) such that in any trading 
(considered as a scattering) the agent saves a fraction A^ of its money nii(t) 
at that time and offers the rest (1 — X^m^t) for random trading. We assume 
each trading to be a two-body (scattering) process. The evolution of money 
in such a trading can be written as: 

rm(t + 1) = Xirriiit) + dj [(1 - X l )m l (t) + (1 - Xj)m 3 {t)] , (1) 

mj(t + 1) = X jmj (t) + (1 - e tf ) [(1 - Xi)mi(t) + (1 - A,)™,(i)] (2) 

where each mi > and is a random fraction (0 < e < 1). In the fixed 
savings market Aj = Xj for all i and j, while in the distributed savings market 
Xi ^ Xj with < A i; Xj < 1. 

3 Numerical observations 

In addition to what have already been reported in Ref. [8, 9, 10] for the model, 
we observe that, for the market with fixed or uniform saving factor A, a fit to 
Gamma distribution [15], 

Q \ 

P(m) ~ to" exp(-m/T), r] = (3) 

1 — A 

is found to be better than a log-normal distribution. However, our observation 
regarding the distribution D(A) of difference A = \Am\ of money between 
any two agents in the market (see Fig. la) suggests a different form: 
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Fig. 1. D(A) in the fixed or uniform savings market, for A = 0.2,0.5,0.8 (right to 
left) and their fitting curves: D(A) ~ exp(— Z\ 1+A /T'); the corresponding P(m) the 
inset. 

P(m) ~ m s exp(-m K /T'); k = 1 + A. (4) 

In fact, we have checked, the steady state (numerical) results for P(m) asymp- 
totically fits even better to (3), rather than to (4). 

With heterogeneous saving propensity of the agents with fractions A dis- 
tributed (quenched) widely (0 < A < 1), where the market settles to a critical 
Pareto distribution Pirn) ~ m~( 1+l/ ) with v ~ I [8], the money difference be- 
haves as D(Am) ~ (zim)~( 1+7 ) with 7 ~ 1. In fact, this behavior is invariant 
even if we set — 1/2 [18]. This can be justified by the earlier numerical 
observation [7, 8] for fixed A market (Aj = A for all i) that in the steady state, 
criticality occurs as A — > 1 where of course the dynamics becomes extremely 
slow. In other words, after the steady state is realized, the third term con- 
taining e = l/2 becomes unimportant for the critical behavior. We therefore 
concentrate on this case in this paper. 

4 Analysis of money difference 

In the process as considered above, the total money (mi + rrij ) of the pair of 
agents i and j remains constant, while the difference Arriij evolves for e = 1/2 
as 

(Am,ij)t+i = aij(Amij) t + Pij{mi + mj) t , (5) 

where = 5 (A.; + Ay) and /3y = 5 (A, — Aj). As such, < a < 1 and 
— i < j3 < 5. The steady state probability distribution D(A) can be written 
as (cf. [18]): 
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D(A) = / dA' D{A') (6(A - (a + 13) A') + S(A - (a - j3)A')) 




where we have used the symmetry of the /3 distribution and the relation 
&ij + fiij — K, and have suppressed labels i 7 j. Here (...) denote average over 
A distribution in the market. Taking now a uniform random distribution of 
the saving factor A, p(A) = 1 for < A < 1, and assuming D(A) ~ 
for large A, we get 

1 = 2 J dX A 7 = 2(1 + 7) _1 , (7) 

giving 7 = 1 . No other value fits the above equation. This also indicates that 
the money distribution P{m) in the market also follows a similar power law 
variation, P(m) ~ to~( 1+1/ ) and v = 7. 



5 Master equation approach 

We also develop a Boltzmann-like master equation for the time development 
of P(m, t), the probability distribution of money in the market [18]. We again 
consider the case = ^ in (1) and (2) and rewrite them as 



J / t+i 



*=(£%)■• "*-5<"» <» 



Collecting the contributions from terms scattering in and subtracting those 
scattering out, we can write the master equation for P(m, t) as 

dP(m t) f f 

\-P(m,t) = { / drrii / drrij P(m i7 t)P(m : j ,t) 5(^+771,+ /j,Jmj-m)}, 

(9) 

which in the steady state gives 

P{m) = ( J drrii J drrij P(mi)P(nij) 5(/i+TO; + ^Jirij — to)). (10) 
Assuming, P(m) ~ m~( 1+l/ ) for m — > 00, we get [18] 

1 = (Oi+r + ( M -n = y J dn+dv-pin+Mn-) + (/x-)"] . (11) 

Considering now the dominant terms (oc x~ r for r > 0, or oc ln(l/x) for r = 0) 
in the x — > limit of the integral / °° m^ u+r "> P{m) exp(—mx)dm, we get from 
eqn. (11), after integrations, 1 = 2/(i/ + 1), giving finally v = \. 
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6 Summary 

We consider the ideal-gas-like trading markets where each agent is identi- 
fied with a gas molecule and each trading as an clastic or money-conserving 
(two-body) collision [7, 8, 9, 10]. Unlike in a gas, we introduce a saving fac- 
tor A for each agents. Our model, without savings (A — 0), obviously yield 
a Gibbs law for the steady-state money distribution. Our numerical results 
for uniform saving factor suggests the equilibrium distribution P(m) to be 
somewhat different from the Gamma distribution reported earlier [15]. 

For widely distributed (quenched) saving factor A, numerical studies 
showed [8, 9, 10] that the steady state income distribution P(m) in the market 
has a power-law tail P(m) ~ m~( 1+l/ ) for large income limit, where v ~ 1.0, 
and this observation has been confirmed in several later numerical studies as 
well [16, 17]. It has been noted from these numerical simulation studies that 
the large income group people usually have larger saving factors [8]. This, in 
fact, compares well with observations in real markets [13, 19]. The time corre- 
lations induced by the random saving factor also has an interesting power-law 
behavior [17]. A master equation for P(m,t), as in (9), for the original case 
(eqns. (1) and (2)) was first formulated for fixed A (Aj same for all i), in [14] 
and solved numerically. Later, a generalized master equation for the same, 
where A is distributed, was formulated and solved in [16] and [18]. We show 
here that our analytic study (see [18] for details) clearly support the power-law 
for P(m) with the exponent value v = 1 universally, as observed numerically 
earlier [8, 9, 10]. 
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